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Objectives 

This  program  focused  on  the  development  of  physics-based  models,  highly  efficient  numerical  approx¬ 
imation  techniques,  parameter  estimation  algorithms,  and  model-based  control  techniques  for  shape 
memory  alloys  and  polymers,  electroactive  polymers,  magnetic  actuators,  and  PZT-based  actuators 
and  sensors  that  provide  unique  transducer  capabilities  for  Air  Force  applications.  In  the  modeling 
and  simulation  component,  energy  and  statistic-based  techniques  were  used  to  characterize  hystere¬ 
sis,  constitutive  nonlinearities,  and  complex  dynamics  in  a  manner  that  facilitates  device  and  control 
design.  The  emphasis  in  the  control  component  focused  on  the  development  of  robust  linear  and 
nonlinear  algorithms  that  permit  real-time  implementation  while  achieving  stringent  stabilization 
and  tracking  specifications. 

Accomplishments 

We  focused  on  shape  memory  alloys,  polymers  and  composites,  electroactive  polymers,  magnetic 
devices,  and  lead  zirconate  titanate  (PZT)-based  devices  due  to  the  unique  actuator,  sensor  and 
structural  capabilities  that  they  provide  for  present  and  future  Air  Force  applications.  Shape  mem¬ 
ory  alloys  (SMA)  are  presently  being  considered  for  vibration  attenuation  in  membrane  mirrors  and 
modular  antennas,  symmetry  breaking  and  flutter  control  in  jet  engines,  chevrons  for  improved  mix¬ 
ing  and  flow  in  jet  engines  to  decrease  noise  and  drag,  and  MEMs  devices  for  flow  control.  Shape 
memory  polymers  (SMP)  provide  unique  deployment,  positioning  and  vibration  attenuation  capabil¬ 
ities  for  payload  booms,  large-scale  mirrors  and  antennas.  Electroactive  polymers  (EAP)  comprised 
of  carbon  nanotube  infused  polyimides  are  being  considered  for  thin-film  membrane  mirrors,  gos¬ 
samer  antennas  and  micro-air  vehicle  (MAV)  skins  and  frames.  Finally,  PZT-based  devices,  such  as 
macro-fiber  composites  (MFC),  are  being  investigated  within  the  context  of  deployment,  positioning 
and  vibration  attenuation  of  large  aerospace  structures  including  thin-film  membrane  mirrors  and 
gossamer  antennas.  However,  all  of  these  compounds  exhibit  highly  nonlinear  and  hysteretic  dynam¬ 
ics  that  must  be  incorporated  in  models  and  model-based  control  designs  before  their  potential  can 
be  fully  realized  in  Air  Force  applications. 

Development  of  the  Homogenized  Energy  Model  (HEM)  Framework  for  Ferroic  Compounds 

Through  previous  AFOSR  support,  we  developed  an  initial  framework  for  characterizing  hystere¬ 
sis  and  constitutive  nonlinearities  in  ferroelectric,  ferromagnetic  and  ferroelastic  materials  which  we 
collectively  term  ferroic  compounds.  Representative  transducer  materials  from  these  material  classes 
are  PZT,  Terfenol-D  and  shape  shape  memory  alloys  (SMA)  such  as  Nitinol.  During  the  current 
program,  this  framework  was  significantly  expanded  to  provide  a  unified  framework  for  characteriz¬ 
ing  the  intrinic  coupling  between  strains,  polarization  and  magnetization.  This  provides  a  natural 
basis  for  subsequent  construction  of  constitutive  relations  for  devices  such  as  beams,  plates,  rods, 
membranes  and  shells  that  incorporate  PZT,  magnetic  or  shape  memory  alloy  actuators  and/or  sen¬ 
sors.  We  summarize  here  the  development  from  [9]  for  piezoelectric  material  to  illustrate  the  general 
framework. 

At  the  lattice  level,  we  consider  the  Helmholtz  and  Gibbs  energy  densities 

^(-P,  e)  =  2??a(-P  -  -Pr)2  +  2C“(e  ~~  £r)2  ha(P  -  Pr){£  ~  £r)  (1) 

and 

Ga{E,  cr;  P,e)  =  ifa(P,£)  —  EP  —  as  (2) 

where  we  indicate  ±180°  and  90°  orientations  by  a  =  ±,90.  Here  P,e,a,  and  E  respectively 
denote  the  polarization,  strain,  stress  and  field  whereas  Pr,  £r,t]£ ,  cp  and  h  denote  the  remanence 
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polarization,  remanence  strain,  inverse  susceptibility  at  constant  strain,  elastic  stiffness  at  constant 
polarization,  and  piezoelectric  constant. 

For  a  fixed  applied  field  E  and  stress  a,  the  conditions  =  0  and  (p-  =  0  yield  the  relations 


K  =  Pr  +  xlE  +  dQ(j 


+  daE  sa<7 


where 


Xa  = 


'la  ILa 


dn  — 


hi  - 


K _  E  =  ffa 

CP  ri£  5  a  71£  —  h 1,^ 

^ a'la  ° a'la  ,0a 


(3) 


(4) 


are  the  ferroelectric  susceptibility  at  constant  stress,  the  piezoelectric  constant,  and  elastic  compliance 
at  constant  field.  The  minimum  of  the  Gibbs  energy  in  each  a-well  can  then  be  expressed  as 


Gam(E,  a)  =  ~\x*E2  -  ^a2  -  daEo  -  EP%  -  ae% 


(5) 


The  dipole  fractions  x+,x_  and  xgo  associated  with  positively,  negatively,  and  90°  degree  dipoles 
evolve  according  to  the  differential  equations 

X-  =  -P-9QX-  +P90-X90 

X90  =  P-90X-  -  (P90-  +  P90+)X90  +  P+90X+  (6) 

x+  =  P90+X90  —  P+9QX+ 


where 


Pap{E,o) 


1  e-A  GZhE,a)V/kT 
Ta/3 


(7) 


quantifies  the  likelihood  of  transitioning  from  an  a-well  to  a  /3-well.  Note  that  (6)  can  be  simplified 
using  the  relation  x+  +  x_  +  xgo  =  1.  The  activation  energy  is  specified  by  the  relation 


f  AGo(l  -  Faf){E,  a)/Fc )2  ,  -Fc  <  Fap(E,  a)  <  Fc 
A Gaa0(E,a;Fc)  =  l  Fc  ,  Fa/3(E ,  a)  <  - Fc  (8) 

[  0  ,  Fa/3(E,  a)  >  Fc. 


Here 

Fap(E,  a)  =  Gam(E,  a)  -  Gpm(E,  a)  (9) 

is  the  thermodynamic  driving  force.  Note  that  A  Go  =  \FC  is  the  energy  barrier  at  zero  driving  force. 
The  polarization  and  strain  kernels  are  given  by 


p=  £ 

ct=±,90 


r  Pa 

Jja±  m  ) 


*  =  £ 

a=±,90 


T  F ^ 
'Xja(~'rn  * 


(10) 


Based  on  (3),  these  relations  can  be  expressed  as 

P(E,  cr)  =  d(E ,  a)a  +  xaE  +  Pirr(E,  a) 
e(E,  a)  =  sEcr  +  d(E ,  a)E  +  £irr(E,  a) 


(11) 


where  __ 

d(E,a)  =  Y  daxa(E,  a) 

a=±,90 

Pirr(E,a)=  Y  PrX<*(E,<j) 

a=±,  90 

£irr{E,a)=  Y  £Rxa{E,  a). 

a=±,90 


(12) 
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In  the  final  step  in  the  development,  macroscopic  models 


roc  roc  _ 

P(E(t),a(ty,x°+)=  /  PiE^  +  E^aityFME^FjdEjdFc 

JO  J — oo 
roc  roc 

e{E(t),a{t)-,x\)=  I  /  e(E(t)  +  EI,a(t)]Fc)iyI(EI)iyc(Fc)dEIdFc 


(13) 


JO  J  —oo 

are  constructed  by  considering  interaction  fields  Ej  and  thermodynamic  driving  forces  Fc  to  be 
manifestations  of  underlying  densities  uj  and  vc.  The  final  models  can  thus  be  expressed  as 


P(E,  cr)  =  d(E,  cr)a  +  x°E  +  Pirr{E,  cr) 
e(E,  cr)  =  sEcr  +  d(E,  a)E  +  £irr(E,  a) 


(14) 


where 


d(E,  a)  = 


'0  J —oo 
roc  roc 


d(Ee ;  Fc)uj(Ej)uc(Fc)dEIdFc 


P }rr(E ,  cr)  — 


'0  J  —oo 
roc  roc 


Pirr{Ee;  Fc)uj(Ej)i'c(Fc)dEjdFc 


(15) 


roc  roc 

£irr(E,a)=  /  /  £irr(Ee;Fc)uI(EI)i'c(Fc)dEIdFc. 

J  0  J —oo 

Details  regarding  various  choices  for  the  densities  are  provided  in  [9].  For  the  implementation 
and  parameter  estimation  algorithms,  we  employ  the  representations 


I<a 
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VC{FC)  =  Cl  akMFC)  ,  MFc)  =  nr- 

k=ka  v*FcV2tt 


-[ln(Fc)-„c]2/2(^)2 


K_^ 

r 

k=k 


(16) 


MEi)  =  C2  PktPkiEO  ,  <Pk{E{)  =  E^/2(<T/fc) 


where 


ci 


\k=ka  J  \k=k/3  J 


(17) 


ensure  integration  to  unity.  During  model  calibration,  the  parameters  {ak,Pk}  are  determined 
through  a  least  squares  fit  to  the  data. 

To  obtain  an  explicit  constitutive  equation,  we  invert  the  strain  relation  (14)  to  obtain 


cr (E,  e)  =  cE£  —  e(E,  <jq )E  —  cE£irr(E,  <to) 


(18) 


where  ^  ^ 

(  ,  e(E,  do)  ^ d(E,  ctq) . 

The  constitutive  relation  (18)  can  be  used  to  construct  rod,  beam  and  shell  models  for  distributed 
structures.  To  illustrate,  consider  a  rod  having  density  p  and  cross-sectional  area  A.  Longitudinal 
displacements  are  denoted  by  u(t,x).  Force  balancing  yields 

d2u  ON 
P  dt 2  dx 

where  the  force  resultant  is  N  =  a  A.  The  stress  is  given  by  (18). 
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Figure  1:  (a)  Aerospace  application  employing  macro-fiber  composites  (MFC).  Model  fit  to  800  VDC 
experimental  data  and  model  predictions  for  100,  200  and  400  VDC  data  in  the  (b)  phase  space  and 
(c)  time  domain,  (d)  Histogram  for  the  bootstrapped  estimate  for  xa ■ 


For  general  elastic  structures  with  coupling  between  the  electric  and  structural  components,  one 
would  employ  the  relations 

V  •  D  =  0  ,  D  =  eo  E  +  P 

pu  =  V  ■  a  (19) 

V  x  E  =  0  ,  E  =  —  \7<p 

with  the  vector  version  of  (18). 

Details  regarding  the  use  of  (18)  to  construct  models  for  PZT-based  macro-fiber  composite  (MFC) 
actuators  operating  in  highly  hysteretic  and  nonlinear  regimes  can  be  found  in  [11].  Representative 
model  fits  are  illustrated  in  Figure  1. 

Data-Driven  Techniques  for  Parameter  Estimation 

To  construct  the  density  basis  functions  ipk(Ej )  and  cj>c(Fc),  it  is  necessary  to  determine  values 
for  the  driving  force  mean  pc  and  variance  a as  well  as  the  interaction  field  variance  a'j.  We  denote 
this  set  of  parameters  by 

P=  [Mc^cV/]-  (20) 

The  remaining  parameters  are  denoted  by 

P  =  [PR,£wER,x‘T,d+,sE,'y,T90,ak,pk}.  (21) 
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For  the  polarization-strain  model,  this  comprises  the  set  that  is  optimized  during  model  calibration. 
Structural  models  require  additional  elastic  and  damping  parameters. 

A  critical  step  when  calibrating  models  is  the  determination  of  initial  parameter  values  that 
facilitate  the  determination  of  optimized  parameter  values  through  least  squares  fits  to  data.  A 
significant  aspect  of  the  program  focused  on  the  development  of  data-driven  techniques  to  determine 
initial  parameter  estimates  based  on  easily-measured  properties  of  polarization,  strain,  or  magneti¬ 
zation  data.  Initial  results  are  reported  in  [1,4, 10]  and  this  comprises  an  active  component  of  the 
present  program. 

Statistical  Model  Development,  Parameter  Estimation,  and  Model  Uncertainty  Quantification 

During  the  program,  we  investigated  the  development  of  statistical  emulator  construction  for 
complex  smart  systems  [18, 19]  and  initiated  the  investigation  of  statistical  techniques  for  parameter 
estimation  and  model  uncertainty  quantification.  It  is  shown  in  [11]  that  standard  use  of  asymptotic 
theory  to  compute  a  covariance  matrix  V  =  a2(xTx)_1i  where  x  is  the  sensitivity  matrix,  is  not 
feasible  since  the  Fisher  information  matrix  xTX  is  singular.  Instead  we  apply  bootstrapping  methods 
to  determine  parameter  densities  and  confidence  intervals.  As  illustrated  in  Figure  1(d),  this  permits 
the  determination  of  non-Gaussian  densities  that  can  be  used  for  subsequent  model  uncertainty 
quantification. 

Model-Based  Control  Design:  Nonlinear  Inverse  Compensators 

The  goal  of  the  control  component  of  the  investigation  was  to  develop  theory  and  algorithms 
that  permit  real-time  tracking  and  vibration  attenuation  using  smart  material  actuators  operating 
in  highly  nonlinear  and  hysteretic  regimes.  For  the  applications  under  consideration,  this  can  result 
in  micron-level  (or  smaller)  tolerances  and  quasistatic  up  to  kilohertz  operating  frequencies. 

As  depicted  in  Figure  2,  there  are  essentially  two  control  strategies  for  hysteretic  systems  using 
nonlinear  models.  The  first  is  to  use  the  characterization  framework  to  construct  an  approximate 
nonlinear  inverse  that  linearizes  the  actuator  response  in  the  manner  depicted  in  Figure  2(a).  Linear 
control  algorithms  are  then  used  to  achieve  control  objectives.  The  second  strategy  is  to  construct 
nonlinear  control  designs  which  yield  input  signals  that  directly  incorporate  actuator  nonlinearities 
as  depicted  in  Figure  2(b).  The  latter  strategy  is  advantageous  when  tracking  trajectories  are  known 
a  priori  whereas  the  inverse  compensators  will  be  required  when  attenuating  unmodeled  or  stochastic 
uncertainties. 

During  the  grant,  we  focused  significant  effort  on  the  development  and  implementation  of  C\ 
adaptive  control  designs  for  systems  employing  PZT  or  magnetic  actuators  and  sensors.  As  detailed 
in  [5-7],  this  development  employs  the  following  steps.  Based  on  the  physical  nature  of  the  hysteresis 
and  constitutive  nonlinearities,  we  first  assume  that  inversion  errors  cr[t)  are  bounded.  Projection- 
based  adaptation  laws  are  then  constructed  using  the  error  x(t)  =  x(t)  —  x(t)  between  the  states 


(a) 


(b) 


Figure  2:  (a)  Linear  control  design  employing  a  nonlinear,  model-based  inverse  compensator  and 
(b)  nonlinear  model-based  control  design. 
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L1  Control  with  Time  Compensation:  r(t)  (blue/solid  line),  x^t)  (red/dashed  line) 


Figure  3:  Performance  of  the  C\  adaptive  control  law  for  tracking  a  reference  signal  that  includes  a 
ramp,  hold,  and  multiple-frequency  component  r(t)  =  2  cos (t)  +  10cos(7rf/5). 


or  outputs  of  the  physical  and  predicted  systems.  Finally,  we  incorporate  low-pass  filters  to  atten¬ 
uate  high  frequency  noise.  To  accommodate  constant,  filter-induced  time  delays,  we  employ  time 
compensation  as  a  final  step  in  the  design. 

Simulation  results  illustrating  the  performance  of  the  method  for  tracking  a  reference  signal  that 
includes  a  ramp,  hold,  and  multiple-frequency  component  r(t)  =  2  cos (t)  +  10cos(7rt/5),  using  a 
hysteretic  PZT-based  transducer,  are  illustrated  in  Figure  3.  It  is  noted  that  errors  are  maintained 
below  a  level  of  1  /jrn  with  the  maximum  errors  occurring  at  points  where  the  derivative  in  the 
tracking  signal  is  discontinuous. 

Details  regarding  additional  inverse  compensator  and  adaptive  control  frameworks  investigated 
during  this  program  can  be  found  in  [2,8,12,22]. 


Model-Based,  Control  Design:  Nonlinear  Perturbation  Control 

The  second  main  thrust  of  the  control  program  has  focused  on  the  formulation,  numerical  imple¬ 
mentation,  and  experimental  implementation  of  nonlinear  perturbation  control  designs  for  systems 
employing  PZT  and  magnetic  actuators  operating  in  hysteretic  regimes.  As  detailed  in  [15-17],  these 
designs  are  developed  in  two  steps:  (i)  construction  of  an  optimal  nonlinear  open  loop  control  signal, 
and  (ii)  construction  of  a  perturbation  feedback  component. 

To  illustrate,  let  r  denote  a  reference  signal  to  be  tracked,  y{t)  =  Cx(t)  denote  observations  and 
let  e(t)  =  Cx(t)  —  r(t)  designate  the  error.  The  augmented  penalty  functional  is  taken  to  be 


pf  r 


J  =  ^[Cx(tf)  -  r(tf)]TP[Cx(tf)  -  r(tf )]  +  T  H-  \T(t)x(t.) 
2  Jt0  1 

where  A  denotes  the  adjoint  variable,  the  Hamiltonian  is 

H  =  ^  [e(t)T Qe{t)  +  uT(t)Ru{t)\  +  XT  [ Ax(t )  +  [H(u)](t)] , 


dt 


and  Q,  R  respectively  penalize  large  errors  and  control  inputs. 


(22) 


(23) 
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Enforcement  of  necessary  conditions  to  minimize  (22)  yields  the  control  input  relation 


(24) 

along  with  the  two-point  boundary  value  problem 

■4 -T 

ST 

II 

•+o 

(25) 

where  z  =  [x,  A]T  and 

F(t,z)  = 

~  Ax(t)  +  [B(u)](t) 

_  -ATX{t)  -  CTQCx(t)  +  CTQr(t)  _ 

(26) 

To  approximate  the  solution  to  (25),  we  employ  a  finite  difference  discretization  defined  on  the 
grid  tj  =  j  At,  where  At  =  ^  and  j  =  0,-  •  ■ ,  N.  Letting  zj  ~  z(tj),  this  yields  the  discrete  system 

^  \zi+ 1  zj\  ~  2  [E(tj,  Zj)  +  F(tj+ 1,  Zj+ 1)] 

E0z0  =  [x0,0]T 

EfzN  =  [0,  —CTPr(tf)]T. 

(27) 

The  solution  of  (27)  can  be  expressed  as  the  problem  of  finding  Zh  =  [zq, 

•  •  • ,  zn\  which  solves 

Hzh)  =  o. 

(28) 

A  quasi-Newton  iteration  of  the  form 

A+l  —  A  i  ck 
zh  ~  zh  +  Sft) 

where  ^  solves 

E'(z^  =  ~E(4), 

(29) 

(30) 

is  then  used  to  approximate  the  solution  to  (28).  The  Jacobian  has  the  form 

“  So  R0 


E\zh)  = 


Si  Ri 


where 


Si  =  -  — 
At 


E0 

I  0 
0  / 


Sn-i  Rn-i 
Ef 


(31) 


A 

-CtQC 


d\B[u*] 

-At 


(32) 


The  representation  for  R{  is  similar. 

It  is  shown  in  [15, 16]  that  an  analytic  LU  decomposition  can  be  determined  for  J-( z fa).  This 
significantly  reduces  memory  requirements  and  is  fundamental  for  efficient  solution. 

The  open  loop  control  provided  by  (24)  is  not  robust  with  regard  to  model  or  measurement 
uncertainties.  To  provide  robustness,  we  have  investigated  two  perturbation  techniques  to  construct 
an  outer  feedback  loop:  PI  and  narrowband  feedback.  The  former  has  the  form 


up i{t)  =  —kpe(t)  —  ki  /  e(s)ds 

Jo 
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Figure  4:  (a)  Reference  trajectory  and  experimental  tracking  performance  obtained  with  PI  control, 
optimal  open  loop  control  and  perturbation  control  at  1000  Hz.  (b)  Hysteretic  and  nonlinear  input- 
output  behavior  of  the  magnetic  actuator. 


where  e  denotes  measured  errors,  and  provides  the  advantage  that  once  the  open  loop  control  has 
been  computed,  it  is  just  as  efficient  to  implement  as  a  straight  PI  design.  We  summarize  in  Figure  4 
representative  results  demonstrating  the  performance  of  this  perturbation  design  for  a  magnetic 
transducer  when  experimentally  implemented  at  1  kHz.  It  is  observed  that  it  provides  a  significant 
improvement  over  PI  control  at  this  frequency  while  operating  in  a  highly  hysteretic  and  nonlinear 
regime.  Further  details  are  provided  in  [15,16]. 

The  second  perturbation  design  employs  a  control  law  of  the  form 

V,(t)  =  U*(t)  +  UNB(t)  +  ui(t)  (33) 

where  ui(t)  provides  integral  feedback  and 

UNB(t)  =  ~[Kf  K][xf  ;  e] 

is  a  narrowband  feedback  term  constructed  by  appending  the  filter  equations 

-jj-  =  AfXf(t )  +  BCx(t ) 

to  the  state.  The  frequency  state  matrix  has  components 

1  0 

where  LOi  is  a  frequency  being  targeted  and  is  the  associated  damping  coefficient;  see  [20,21]  for 
details  regarding  the  manner  in  which  hysteretic  actuators  excite  higher  harmonic  responses. 

The  performance  of  this  design  for  a  magnetostrictive  milling  device  employed  at  Etrerna  Prod¬ 
ucts,  Inc.,  is  illustrated  in  Figure  5.  Figure  5(a)  demonstrates  the  model  fit  for  various  minor  loops 
with  the  largest  having  displacements  on  the  order  of  460  /xm.  The  experimental  tracking  errors 
achieved  when  feeding  back  on  various  frequency  ranges  is  illustrated  in  Figure  5(b).  Following  an 
initial  transient  phase,  the  perturbation  feedback  design  (33)  provides  tracking  errors  on  the  order 
of  10  fim,  thus  demonstrating  the  capability  the  design  for  high-accuracy  tracking  in  high-frequency 
regimes.  Further  details  are  provided  in  [17].  This  also  represents  one  aspect  of  the  technology 
transfer  associated  with  the  program. 
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Figure  5:  (a)  Displacement  data  and  fit  of  the  homogenized  energy  model,  (b)  Tracking  performance 
provided  by  the  perturbation  design  (33)  when  feeding  back  on  100  Hz,  100-300  Hz,  and  100-700  Hz 
higher  harmonic  excitations. 
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Transitions/Interactions 

Transitions 

1.  Magnetostrictive  Actuators  -  Etrema:  The  nonlinear  magnetic  constitutive  models  and  pertur¬ 
bation  control  techniques  developed  in  part  through  AFOSR  support  were  extended,  in  collab¬ 
oration  with  scientists  at  Etrema  Products,  to  optimize  the  performance  of  Active  Machining 
Systems  (AMS)  used  to  mill  products  such  as  piston  heads  [17].  The  goal  is  to  significantly 
increase  milling  speeds  (e.g.,  up  to  20,000  cycles  per  second)  while  maintaining  micron-level 
accuracy.  The  Etrema  and  AFOSR  investigations  were  mutually  complementary  and  each  pro¬ 
vides  substantial  technology  transfer  to  the  other.  Point  of  Contact:  Rick  Zrostlik,  Etrema 
Products  Inc.,  Ames,  IA,  515-296-8030. 

2.  Shape  Memory  Polymers  -  AFRL:  We  investigated  the  development  of  models,  simulation 
techniques,  and  control  designs  for  shape  memory  polymers  presently  being  considered  by 
AFRL  researchers  for  deployment  of  gossamer  structures.  The  shape  memory  polymers  provide 
the  work  density  requirements  necessary  to  deploy  aerospace  structures  such  as  radar  antennas 
and  large  optical  mirrors  but  exhibit  viscoelastic  and  thermodynamic  attributes  that  are  not 
accommodated  by  present  models  or  control  techniques.  Point  of  contact:  Thomas  Murphey, 
Kirtland  AFB,  Albuquerque,  NM,  505-846-9969. 
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Conference,  Colloquia  and  Workshop  Presentations 

•  Eighteenth  International  Symposium  on  Mathematical  Theory  of  Networks  and  Systems  (MTNS 
2008),  Blacksburg,  VA,  July  29,  2008. 

•  Colloquium,  Department  of  Electrical  and  Computer  Engineering,  University  of  Wyoming, 
Laramie,  WY,  September  12,  2008. 

•  Applied  Math  Seminar,  Duke  University,  Durham,  NC,  September  22,  2008. 

•  SPIE’s  16th  Symposium  on  Smart  Structures  and  Materials,  San  Diego,  CA,  March  9,  2009. 

•  2009  U.S.  Navy  Workshop  on  Acoustic  Transduction  Materials  and  Devices,  Penn  State  Uni¬ 
versity,  State  College,  PA,  May  13,  2009. 

•  SIAM  Conference  on  Control  and  Its  Applications,  Denver,  CO,  July  7,  2010. 

•  SPIE’s  17th  Symposium  on  Smart  Structures  and  Materials,  San  Diego,  CA,  March  8,  11  2010 
(3  presentations). 

•  Colloquium,  Department  of  Mathematics,  University  of  Waterloo,  Waterloo,  Ontario,  June  15, 

2011. 

•  Colloquium,  Applied  Mathematics  Department,  University  of  Colorado  Boulder,  August  27. 

2011. 

•  SPIE’s  18th  Symposium  on  Smart  Structures  and  Materials,  San  Diego,  CA,  March  7-8k,  2011 
(McMahan,  Hu). 

•  Symposium  on  Hysteresis  Modelling  and  Micronragnetics  (HMM  2011),  Levico,  Italy,  May  11, 

2011. 

•  Invited  Presentation,  7th  Workshop  on  Control  of  Distributed  Parameter  Systems  (CDPS 
2011),  Wuppertal,  Germany,  July  20,  2011. 

•  Invited  Presentation,  SIAM  Conference  on  Control  and  Its  Applications,  Baltimore  MD,  July 
25,  2011. 


Publictions  Publications  resulting  from  work  supported  by  this  grant  are  listed  as  references. 
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